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ABSTRACT 

The outer (>30 AU) regions of the dusty circumstellar disk orbiting 
the ~2-5 Myr-old, actively accreting solar analog LkGa 15 are known 
to be chemically rich, and the inner disk may host a young protoplanet 
within its central cavity. To obtain a complete census of the bright¬ 
est molecular line emission emanating from the LkGa 15 disk over the 
210-270 GHz (1.4 - 1.1 mm) range, we have conducted an unbiased radio 
spectroscopic survey with the Institute de Radioastronomie Millimetrique 
(IRAM) 30 meter telescope. The survey demonstrates that, in this spec¬ 
tral region, the most readily detectable lines are those of GO and its iso- 
topologues ^^GO and G^®0, as well as HGO+, HGN, GN, G 2 H, GS, and 
H 2 GO. All of these species had been previously detected in the LkGa 
15 disk; however, the present survey includes the hrst complete coverage 
of the GN (2-1) and G 2 H (3-2) hyperhne complexes. Modeling of these 
emission complexes indicates that the GN and G 2 H either reside in the 
coldest regions of the disk or are subthermally excited, and that their 
abundances are enhanced relative to molecular clouds and young stel¬ 
lar object environments. These results highlight the value of unbiased 
single-dish line surveys in guiding future high resolution interferometric 
imaging of disks. 
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1. Introduction 


To understand the origins of our own solar system, it is necessary to investigate 
the protoplanetary disks of young, solar-mass stars (T Tauri stars) that represent 
early solar analogs. The LkCa 15 star/disk system, which has been extensively 
studied at X-ray, optical, infrared, and (sub)millimeter wavelengths, represents a 


prime example of such an early solar analog (Thi et al. 2004 Oberg et al. 


2010 


Andrews et ah 1 120111 llsella et al. 1 120121 ISkinner & Giidell 120131 IThalmann et al 


2014). The LkCa 15 system is associated with the Taurus-Auriga star-forming cloud 


complex, at a distance of about 140 pc (|van den Ancker et al. 1998|), and is an 
archetypal 


“transition disk” system (see, e.g., Williams & Cieza 2011[ ). The K5 
{Lhoi = 0.74Lq) central star is estimated to have an age of 2-5 Myr and a mass 


of M = 1.0 Mq (Simon et al. 2000; Kenyon & Hartmann 1995). The LkCa 15 
star/disk system is characterized by a mass accretion rate of about 3 x 10“® M(t 


© 


yr ^ (jHartmann et al. 

1998 

), a partially dust-depleted inner region of ~50 AU in 

radius ( 

Andrews et al. 

2011 

Espaillat et al.||2007 

Pietu et al.||2006 

), and an infrared 


excess over the stellar photosphere that can be explained by the presence of hot dust 


within a few Alls from the central star (Espaillat et al. 2007 Mulders et al. 2010) 


The disk, which displays detectable molecular gas in Keplerian rotation out to ~900 
AU (Pietu et al. 2007), is viewed at intermediate inclination {i ~ 52°; Dutrey et al. 


2011 ). 


The ~50 AU central cavity constitutes one of the largest inner submm “holes” 
observed among transition disks ( Andrews et al.||2011 ). The discrepancy in the outer 
disk radius as inferred from interferometry of emission from dust and gas (150 AU 
vs. 900 AU, respectively; Pietu et al.||2007 Isella et ^|2012 ) also suggests the density 
of mm-size dust grains drops precipitously in the outer disk, beyond ~150 AU. This 
radially stratihed dust disk structure may reflect the effects of recent or ongoing 


planet formation. Indeed, Kraus & Ireland (2012) reported the potential discovery 
of a 6 Mj, ~ 1 Myr old protoplanet well within the submm continuum cavity (i.e., at 
a distance of 16 AU from the central T Tauri star). It is unclear how a cavity as large 
as that observed within the LkCa 15 disk may be formed, but it is believed that it 
must be due to more than just dynamical interactions with a single planet, suggesting 
the presence of multiple planets and/or the need for additional clearing mechanisms 


(e.g., grain growth, photoevaporation; Zhu et al. 2011; Dodson-Robinson & Salyk 
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2011 ). 


Given its status as the first accreting T Tauri star that has been identihed as 
potentially harboring a young protoplanet, the LkCa 15 system serves as a unique 
laboratory for constraining physical conditions within a planet-forming disk — one in 
which planet formation is likely at an advanced stage. Indeed, among protoplanetary 


disks, the LkCa 15 disk is especially chemically rich; it has 


in several molecular transitions (Thi et ah 


Pietu et al. 


Dreviously been detected 
Oberg et al. 1|2010| ). 


2007 


Furthermore, Skinner & Giidel (2013) report the detection of LkGa 15 as a bright 
X-ray source with Ghandra (intrinsic X-ray luminosity of log ~ 30.4 ergs s“^). As 


discussed in Kastner et al. (2014) and references therein, such strong X-ray emission 


(as well as extreme untraviolet emission) from the central star likely has important 
consequences for disk heating, chemistry, accretion, and mass-loss — and, hence, for 
planet formation. 

We have demonstrated that unbiased line surveys of evolved, chemically rich, 
protoplanetary disks, such as that orbiting LkGa 15, are necessary to establish the full 
inventory of readily detectable molecular species in disks ( Kastner et al.||2014 ). Such 
observations are important both in providing constraints on the latest generations 
of disk models that incorporate stellar irradiation and gas-grain reactions, and in 
serving as guidance for subsequent interferometric imaging of molecular emission 


from protoplanetary disks (Kastner et al. 2015). With this as motivation, we have 


carried out such an unbiased radio molecular line survey of the LkGa 15 disk with 
the 30 m IRAM radio telescope. 


2. Observations and Data Reduction 

The line survey of LkGa 15 reported here was carried out over the frequency 
range ~ 210 — 270 GHz with the IRAM 30 m telescope and its EMIR receivers in 
March 2012 during good weather conditions (~ 2.0 mm of precipitable water vapor). 
The EMIR receiver temperatures were fairly homogeneous over the frequency range 
surveyed (see Figure [^, with values in the range 80-100 K. The standard three- 
phase calibration procedure at the IRAM 30 m was used to obtain the zenith opacity 
every 10-15 min, with measured values between 0.1 and 0.3. The EMIR receivers 
were used along with the Fourier Transform Spectrometers in their 200 kHz (0.25 
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km s“^) spectral resolution mode, providing ~ 16 GHz instantaneous bandwidth 
that was split into two chunks of ~ 8 GHz each (see Figure [^. To ensure we 
obtained flat baselines and optimal atmospheric subtraction, the observations were 
performed in the wobbler switch mode, with a 70” throw. The pointing was checked 
approximately every 1-1.5 hours on a nearby continuum source. The focus was 
checked on a strong source every 2-2.5 hours. The frequency setups were chosen 
to ensure a factor two redundancy, except at the edges. However, bad weather 
conditions reduced the observing time, and the hnal frequency coverage was reviewed 
accordingly (see Figure]^. Gonsequently, 9 setups have been observed so as to obtain 
a complete coverage from 208.546 to 269.666 GHz. We note that having a large 
instantaneous bandwidth not only minimizes inter-calibration issues between lines, 
but also minimizes pointing issues. 

Data reduction was performed with the GLASS software within the GILDAS 
packag^ Residual bandpass effects were subtracted by applying a hrst order poly¬ 
nomial to each spectrum. We hnd the noise to be almost uniform (see Figure [^, with 
an average noise level of 7 mK (antenna temperature scale) in 0.27 km s“^ channels 
across the full bandwidth. 

The output of the foregoing data acquisition and reduction sequences for LkGa 
15 is a single spectrum that covers the full ~ 210-270 GHz frequency range. The 
results for line intensities and derived quantities described in Section were obtained 
from these spectra, adopting a beam full width at half maximum (FWHM) of 10.7”. 
The line intensities are quoted in main beam temperature scale, using T^b = 
where the appropriate values of ^eff (forward efficiency) and i?e// (main beam ef- 
hciency) were found by linear interpolation of the values that are provided in the 
IRAM 30 m documentation. These line intensities can be converted to Jy km s“^ 
using 3.906^^ Jy/K for S'/T^, the point source sensitivity, with F^jf = 0.88 — 0.94 
and Agff = 0.41 — 0.49 (aperture efficiency), where these represent the values for the 
IRAM 30 m telescope in the 1.3 mm windo'vj^ 


^ http: / / www. iram. fr/IRAMFR/ GILDAS 

^See http://www.iram.es/IRAMES/mainWiki/IramSOmEfficiencies. 
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3. Results and Analysis 
3.1. Molecular Line Inventory 


We used the spectral line identification methods available in the WEEDS exten¬ 
sion of the GILDA^ software tools, along with visual inspection, to compile a list of 
molecular transitions that are readily detectable and measurable in the 210-270 GHz 
IRAM spectrum of LkGa 15. A summary of these molecular transitions is presented 
in Table and Table The spectral regions covering all the lines we detected are 
displayed in Figure]^ A listing of all molecular transitions previously detected and 
searched for in LkGa 15 is presented in Table An inspection of Table shows 
that, in the spectral region surveyed, the most readily detectable lines are those of 
GO and its isotopologues ^^GO and G^®0, as well as HGO+, HGN, GN, G 2 H, GS, 
and H 2 GO. All of these species had been previously detected in the LkGa 15 disk. 
In most cases, however, the detections were widely spaced in time and/or made 
with different telescope/instrument combinations, potentially resulting in significant 
calibration uncertainties where, e.g., measurement of line ratios is concerned. 


Absent such uncertainties we can, e.g., exploit our near-simultaneous coverage 
of GO and its isotopologues ^^GO and G^®0 to estimate the optical depths of the 
observed transitions of GO and ^^GO (§ 3.4). In addition, the survey includes first 


measurements of the full suite of hyperfine transitions of GN A = 2 —)■ 1 and G 2 H 
A = 3 —)■ 2, whose analysis yields estimates of optical depth and determinations of 
excitation temperature for these molecules 


3.5). 


We note that the J = 5 —)■ 4 
transition of GS is detected here for the first time in the LkGa 15 disk. We also 
detect multiple transitions of H 2 GO, one of which (the J = 3i3 —?• 2 i 2 line at 211.211 
GHz) had not been previously detected in this disk. 


3.2. Line Intensity Measurements 


Because the LkGa 15 disk is viewed at intermediate inclination {i ~ 52°; Dutrey 


et ah 2011), it displays molecular line profiles that are double-peaked. Thus, we fit 


^See 


Maret et al. 


2011 
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the observed profile with a Keplerian disk model described in Section 2 of Beckwith 


& Sargent (1993) and used by Kastner et ah (2008, 2010, 2014) to analyze molecular 
emission of protoplanetary disks observed with a single-dish radio telescope. The 
prohle parameters are: the peak line intensity {Tmb), the Keplerian velocity at the 
outer edge of the disk equivalent to half of the peak-to-peak velocity width), the 
disk radial temperature prohle power law component {q, where T(r) a and the 
outer disk density cutoff (pd). When we have obtained a good ht to the observed line 
prohle, as determined by the reduced of th® fit) 'w® can then integrate under the 
model line prohle and obtain the total line intensity (/). We stress that the purpose 
of applying such a (parametric) model is to obtain line intensity measurements, and 
not necessarily to infer disk structure; such inferences are, in any case, inherently 
difficult given only single-dish observations. 


For all lines with high signal-to-noise, all four parameters were allowed to vary 
to hnd the best-htting line prohle, as determined by minimizing the reduced of the 
ht. To ht the low signal-to-noise lines, we adopt the model parameters for Vd, q, and 
Pd from the line and only allow T^b to remain free. To estimate upper limits 

and tentative detections, we hxed all parameters, except Tmb, to the values that best 
ht the line, and then determined the maximum line intensity that could be 

present given the observed noise level. The upper limit (at the one-sigma level) was 
then determined by integrating over the resulting model line prohle. The results 
of these line prohle hts are listed in Table [T] the upper limits listed are at the 3-a 
level. In Figures |4]j^ we display the observed and best-ht model line prohles for the 
transitions of CO isotopologues, HCO’*', HCN, CS, and H 2 CO; in all of these cases, 
all four model prohle parameters were allowed to vary. We also obtain tentative 
detections for DCN and DCO’*', for which 3 of the 4 parameters in the ht were hxed. 
The spectra of these tentative detections were rebinned and are presented in Figure 
For comparison, a summary of previous molecular line intensity measurements for 
LkCa 15 is provided in Table 


3.3. Column Density Estimates 

We estimate source-averaged column densities for the (pure rotational) molecu¬ 
lar transitions detected using the methods described in Goldsmith & hanger (1999) 
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and [Kastner et aL (2014), where we use the results for the integrated line intensities 
(Table [^, along with molecular line reference data from The Cologne Database for 
Molecular Spectroscopy (CDMS)|^ These methods are briefly summarized here. 


In the case of optically thin emission (r ^ 1), the total column density for 
linear molecules - for which the degeneracy of the upper rotational level is given 
by gu = ‘^Ju + 1 and the energy of the lower level is Ei = hBrotJiiJi + 1) ^ can be 
approximated by 


Ntot = 


QrotjTex) 
-^ulQu 1 


^Ei/kTea: PP 

A T fdilutions 


( 1 ) 


where u is the frequency of the transition, W is the integrated line intensity (in K 
cm s“^), Aui is the Einstein coefficient for the transition, T^i = ^ is the temperature 


equivalent of the transition energy, fdilution = 




is the main-beam dilution cor¬ 


rection factor (with 12^6 and the main-beam and source solid angles, respectively), 
and AJ^ = J^iTex) - Ju(Tbg) where J^, = x exp[hu)kT)-i T^g is the temperature 
of the cosmic microwave background radiation. For optically thick emission (r ~ 1), 
a correction factor of must be applied to Equation 


For the simple (diatomic) case where the rotational level is not split by higher- 
order effects (e.g., fine, hyperfine structures) the calculation of the total column 
density is straightforward. The temperature dependent total partition function then 
reduces to the rotational partition function, which can be approximated by 


Z{Tex) ~ Qtot{Tex) ~ QrotiTex) ^ 


1 

3 


kE, 


hB. 


rot 


( 2 ) 


where Tex is the excitation temperature of the gas and Brot is the molecular rotational 
constant (see, e.g., Mangum & Shirley 2015[ ). Equationonly applies to the case of 
linear (and approximately linear) molecules, such that the molecules have only one 
(or only one dominant) moment of inertia, or cases in which the hyperfine structures 
are not resolved. In the (more complex) case of the H 2 CO molecule, we adopt the 
values of the partition function provided by the CDMS database and interpolate to 
determine the value for the partition function at arbitrary temperatures. 


^See http://www.astro.uni-koeln.de/cdms. 




















Using the foregoing methods, we obtain source-averaged column densities Nx 
for a range of excitation temperatures {Tex = 4, 6, 10, 20, 30 K) for the species 
that are assumed to be optically thin. Molecular data (e.g., v, Brot, Aui) were 
obtained from the CDMS database. In Table we list these column densities for the 
foregoing excitation temperatures, along with column density estimates for (2- 
1). Although the (2-1) emission is evidently optically thin (§ 3.4), the column 

density estimates for (2-1) were derived from the column density obtained 

from the integrated intensity, adopting isotopic abundance ratios ~ 68 

(Milam et ah 2005) and ~ 480 (Scott et al. 2006). 


For all molecules except CN and C 2 H (see § 3.5), the values of column density 


(Table were obtained assuming the characterisitic emitting region is a uniform 
disk of angular radius R = 4", based on results from interferometric studies (e.g., 
Pietu et ^[2007 ; Oberg et al.||2010 Chapillon et al.||2012a ). Hence, we have adoptee 
a beam dilution facto]0of fdilution ~ 1-8 for all molecules in the disk, corresponding 
to the half power beam width of the IRAM 30 m telescope at 230 GHz (10.7”). 
This simplihcation introduces modest systematic errors in the values of fdilution and, 
hence, column density at the extremes of the frequency range surveyed (i.e., a ~ 6% 
underestimate at 210 GHz and a ~ 30% overestimate at 260 GHz). 


3.4. Opacities in CO Lines 

We can take advantage of the detection of G^®0 to estimate the optical depths 
of the ^^GO (2-1) and ^^GO (2-1) transitions. The optical depth of ^^GO may be esti¬ 
mated by comparing the expected and measured line ratios for the GO isotopologues. 
We expect 

^ 1 - exp (-ri3co) 

l-exp(-ri3coA)’ ^ ^ 

where R is the measured ^^GO:G^®0 line ratio and X is the ^^GO:G^®0 abundance 
ratio. Equation follows directly from the equation of radiative transfer, if it is 
assumed that the contribution from the background radiation is negligible and the 
excitation temperatures of ^^GO and G^®0 are identical. We adopt a value of X = 7 


Tn 


Kastner et al. 


(2014), we mistakenly quote the inverse of the beam dilution factor. 
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for the abundance ratio, based on the values = 68 and 

= 480 (as determined from measurements of CN and CO isotopologues 


for the local interstellar medium and the Solar System, respectively; [Milam et al 
2005; Scott et al. 2006). So, if the emisson from is optically thin {tuco 

1), we would expect to obtain R ^ 7. The measured ratio for the J = 2 —)■ 1 
transitions for LkCa 15, R ~ 6.0, therefore conhrms that the emission in the 
(2-1) line is optically thin, with an implied optical depth of tisqq ~ 0.38. Similar 
results were obtained by Qi et al. (2003) from observations of LkCa 15 in the J = 
1 —)■ 0 transitions of and C^®0. Then, given the assumptions for the isotopic 
abundances of carbon, we expect ~ Q^tisco- Hence, we estimate an optical 

depth of ~ 26 for (2-1), i.e., the emission in this transition of from the 

LkCa 15 disk is highly optically thick. 


We caution that the foregoing estimates of ^^CO optical depths are potentially 
subject to large uncertainties. These CO opacity estimates assume that there is no 
selective photodissociation or chemical fractionation across the line prohle. Finally 
(and perhaps least signihcantly), the ^^CO line prohle is slightly asymmetric (see 
Figure]^, which suggests that the telescope was mispointed by a few arcseconds for 
that particular observation, and hence, the intensity of the ^^CO line may be a slight 
underestimate. 


3.5. CN and C 2 H Hyperfine Component Analysis: Optical Depths, 
Excitation Temperatures, and Column Densities 


The CN and C 2 H rotational emission lines display hyperhne splitting (HFS) due 


to interactions between electron and nuclear spins (e.g., Ziurys et al. 1982, see Ta¬ 
ble §. We performed measurement and analysis of the LkCa 15 CN iV = 2 —)■ 1 and 
C 2 H A = 3 —)■ 2 spectra — which are displayed in Figurej^and Figure]^ respectively 
— using a new adaptation of the classical HFS htting method not yet implemented 
in the CLASS software. This adaptation uses two Gaussians of the same width to 
ht the double-peaked emission features in our spectrum that result from the inclined 
Keplerian disk. The two-Gaussian HFS code makes the simplifying assumption that 
all the hyperhne transitions within a given rotational transition share the same exci¬ 
tation temperature. Fitting all the hyperhne components simultaneously results in 
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the determination of the opacity, r, and excitation temperature, Tea,, of the rotational 
transition (see Kastner et ah 2014). Figure]^ Figure]^ and Tablesummarize the 
results of this analysis for the IRAM 30m measurements of the CN and C 2 H emission 
from LkCa 15. 

The determination of the total partition functions Qtot for the CN and C 2 H 


radicals must account for the nuclear-spin statistics, as discussed in Kastner et ah 


(2014). The column density is then computed as 

jkTex 


Ntot = 


Svrz/'* 


Tuldu 


Qtot(Tex) 


( 4 ) 


* 1 - e-^^dTex J ' AulQu 

where all parameters refer to the particular hyperhne transition of interest. We 
assume the characteristic emitting region is a uniform disk of angular radius R = 5", 


based on results from interferometric studies (e.g., Chapillon et al. 


2012a 


Oberg 


et al. 2010). We adopt a beam dilution correction factor fdilution ~ 1-1, but the 


results we obtain are fairly insensitive to the assumed source solid angle. Following 
this procedure, we obtain values of T^x ~ 3.9 K and r ~ 2.2 for CN, and T^x ~ 2.9 
and r ~ 18 for C 2 H, leading to estimates ofNtoj(CN) ~ O.SxlO^'^ cm“^ and Ntof(C 2 H) 
~ 2x 10^® cm“^ for the LkCa 15 disk (Table|^. The results for r indicate the emission 
in both lines is optically thick, where the reader is reminded that the resulting values 
of r represent the total opacity over the HFS. Hence, the derived values of Tgx 
imply either that the gas is subthermally excited or that the CN and C 2 H emission 
originates from disk gas at very low temperature. We caution, however, that the 
low Tex for CN and C 2 H rest on the assumed beam hlling factor. In addition, a less 
extended disk in these molecules would imply a higher Tex- 


Discussion 


In Table]^ we list fractional molecular abundances relative to ^^CO, N(X)/N(^^CO), 
as obtained from the calculated source-averaged molecular column densities for LkCa 
15 


3.3). We present values of N(X)/N(^^CO) that we obtained by adopting Tex = 


4, 6, 10, 20, 30 K. We caution that these column density estimates may have large 
systematic uncertainties, given the uncertainties in the assumed source parameters 
and, in particular, the assumed source emitting region solid angles in the various 


molecules (Pietu et al. 


2007 


Oberg et al. 


2010 


Aikawa et al. 


2003 


Chapillon et al. 
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2012a). Nevertheless, our results yield column densities that are overall consistent 


with previous studies of the LkCa 15 disk. This consistency is apparent from Ta¬ 
ble]^ in which we present a comparison of our column densities and column densities 
determined by previous studies. In most cases, our column density estimates agree 
with previous estimates to within a factor ~3. Apart from CN and C 2 H (see below), 
the largest discrepancies appear in the comparison with results for ^^CO, C^®0, and 


HCO"'" published by Qi et ah (2003); the large differences between column density 


results in these cases can be attributed to the constrasting beam sizes of OVRO 
(~ 3”) and IRAM (~ 12”). 

It should be noted that for LkCa 15, we derive column densities for CN and espe¬ 
cially for C 2 H that are relatively high; in the latter case, the inferred source-averaged 
column density is similar to that of ^^CO. Specihcally, we derive N(CN)/N(^^CO) 
and N(C2H)/N(^^C0) values of approximately 0.01 and 3.4, respectively, adopting 
Tea; = 4 K. The former result is somewhat smaller than that quoted for LkCa 15 in 


Kastner et al. (2014), where the difference reflects our rehned treatment of CN optical 


depth and excitation temperature (§ 3.5), as well as different assumptions regard¬ 


ing beam hlling factors. In deriving our column densities, we assumed that all the 
molecular emission is distributed within a smooth disk, but clearly this may not be 
the case. If the CN or C 2 H emission were instead restricted to specihc regions of the 
disk, we could be severely underestimating the column densities of these molecular 
species. Thus, it is imperative to obtain high signal-to-noise interferometric images 
to study the surface brightness distribution for these highly abundant molecules. 

It is possible that high-energy (FUV and/or X-ray) radiation from the central 


star may be enhancing the abundances of CN and C 2 H in T Tauri disks (Kastner 


et al. [1997 , 2008, 2014 Dutrey et al. [1997 ; Thi et al. 2004 Henning et al. 2010). 


Suggestions about the influence of X-rays on disk chemistry were made by Kastner 


et al. 

(2008 

), 

Thi et al. 

2004 

), 

Salter et al. 

(2011 

) to explain enhanced HCN and 

HCO’*' line fluxes, and bj* 

Kastner et al. 

(2014 

) to explain the enhanced abundances 


that were found for CN and C 2 H for TW Hya and V4046 Sgr. 

We also hnd very low values of excitation temperature for CN and C 2 H (Tea, ~ 
3-5 K) in the LkCa 15 disk. These determinations of low may indicate that the 
emission from CN and C 2 H arises from regions that are deep within the disk. Hence, 


as noted by Kastner et al. (2014), it may be that the abundances of CN and C 2 H 
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are being enhanced by X-rays as opposed to extreme FUV photons, since X-rays can 


penetrate deeper into disks (e.g.. 

Skinner & Giidel 

2013 

Gleeves, Adams, & Bergin 

2013 

). Indeed, 

Skinner & Giidelf 

(2013 

) found that LkGa 15 is an X-ray luminous 


source, with a spectral energy distribution that contains contributions from both soft 
and hard components. This combination may increase the abundances of species (like 
C 2 H) that are the dissociation products of complex organic molecules. Alternatively, 
the CN and C 2 H molecules may reside in tenuous regions of the disk where collision 
probabilities are low, such that the molecules are subthermally excited. 


Our estimates of Ntoi(CN) and Ntof(C 2 H) for the LkCa 15 disk (Table are 


similar to those found for the TW Hya and V4046 Sgr disks (Kastner et ah 2014 


their Table 5). Hence, it appears that large relative abundances of C 2 H (as well as 
CN; see Chapillon et ah |2012a ) may be a common feature of evolved protoplanetary 
disks. The excitation temperatures that we hnd for CN and C 2 H {T^x K) 

in the LkCa 15 disk are similar to, though even lower than, those determined by 


Kastner et al. (2014) for TW Hya and V4046 Sgr; futhermore, we have inferred 


relatively large values for the optical depths of CN and C 2 H for all three of these 
disks. Subsequently, Submillimeter Array imaging of TW Hya revealed that C 2 H 
emission emenates from a ring-like structure, with an inner radius (~ 45 AU) that 
appears to he between the CO “snow line” (as traced by the ~ 30 AU-diameter 
hole in N 2 H^) and the outer boundary of large-grain emission at ~ 60 AU (Kastner 
2015[). Given the similarity of the TW Hya and LkCa 15 disks in terms of 


et al. 


the excitation and optical depth of C 2 H, it is likely that the C 2 H (and possibly CN) 
emission from LkCa 15 also originates from a ring-like structure. 


We present fractional abundances with respect to ^^CO, N(X)/N(^^CO) in Ta¬ 
ble where we cite the results of Thi et al. 2004 as a cross-check of our results. 
This table illustrates that CN and C 2 H are enhanced in evolved T Tauri star disks 
relative to other astronomical objects, especially with respect to YSOs and even the 
expanding, C-rich envelope of the classical carbon star IRC-l-10216, in which the 
abundances of carbon-bearing molecules are strongly enhanced. Thus, our line sur¬ 
vey of LkCa 15 reinforces previous assertions (e.g., Kastner et al.||2014 and references 


therein) that there must be some mechanism that is enhancing the abundances of 
CN and C 2 H in protoplanetary disks. 
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5. Summary 


We have used the Institute de Radioastronomie Millimetrique (IRAM) 30 meter 
telescope and the EMIR fast fourier transform spectrometers to conduct a compre¬ 
hensive mm-wave emission line survey of the circumstellar disk orbiting the nearby, 
pre-main sequence (T Tauri) star LkCa 15 over the 210-267 GHz (1.4 - 1.1 mm) 
range. We hnd that lines of ^^CO, HCO’*', HCN, ^^CO, CN, C 2 H, CS, H 2 CO, and 
C^®0 constitute the strongest molecular emission from the LkCa 15 disk in the spec¬ 
tral region surveyed. The J = 5 —)■ 4 transition of CS and the (J = 3 i 3 —t 212 ) 
transition of H 2 CO are detected for the first time in the LkCa 15 disk. We use 
simultaneous measurements of CO and its isotopologues (i.e., ^^CO and C^®0) to 
estimate the optical depths of the transitions of CO and ^^CO observed; we hnd that 
^^CO emission is optically thin, and emission from ^^CO is highly optically thick. 
The survey also includes hrst measurements of the full suite of hyperhne transitions 
of CN N = 2 ^ 1 and C 2 H = 3 —)■ 2. Modeling of these CN and C 2 H hyper¬ 
hne complexes indicates that the emission from both species is optically thick and 
suggests either that the emission from these molecules originates from either very 
cold (<10 K) regions of the disk, or that the molecules reside in less dense upper 
atmosphere regions of the disk and are subthermally excited. We also hnd that the 
column densities for CN and C 2 H are comparable to those of isotopologues of CO, 
implying that these molecules are very abundant in the disk of LkCa 15. 


As with our molecular line surveys of TW Hya and V4046 Sgr (Kastner et ah 


2014), the results of our molecular line survey of LkCa 15 stress the need for addi¬ 


tional molecular line surveys of protoplanetary disks that cover a range of pre-MS 
evolutionary states; interferometric imaging of protoplanetary disks in transitions 
of CN and C 2 H, as well as other potential molecular tracers of UV and/or X-ray 
radiation; and detailed modeling of T Tauri star-disk systems that includes the ef¬ 
fects of high-energy radiation on disk chemistry and is constrained by observations 
of pre-MS UV and X-ray helds. In general, we note that interferometric spectral 
surveys will become very important in the upcoming years. Unbiased line surveys 
carried out with interferometers would allow for the simultaneous determination of 
the chemical inventory of disks and the emitting regions sizes and morphologies of 
these molecules. 




- 14 - 

This work is based on observations earried out with the IRAM 30m Teleseope. 
IRAM is supported by INSU/CNRS (Franee), MPG (Germany) and IGN (Spain). 
This researeh is supported by National Seience Foundation grant AST-1108950 to 
RIT. The authors thank Gharlie Qi for helpful diseussions. 

REFERENCES 

Aikawa, Y., Momose, M., Thi, W.-F., et al. 2003, PASJ, 55, 11 

Andrews, S. M., Rosenfeld, K. A., Wilner, D. J., & Bremer, M. 2011, ApJ, 742, L5 

Beckwith, S. V. W., & Sargent, A. 1. 1993, ApJ, 402, 280 

Blake, G. A., Sandell, G., van Dishoeck, E. F., et al. 1995, ApJ, 441, 689 

Ghapillon, E., Guilloteau, S., Dutrey, A., Pietu, V., & GuGin, M. 2012a, A&A, 537, 
A60 

Ghapillon, E., Dutrey, A., Guilloteau, S., et al. 2012b, ApJ, 756, 58 

Gleeves, L. I., Adams, F. G., & Bergin, E. A. 2013, ApJ, 772, 5 

Dodson-Robinson, S. E., & Salyk, G. 2011, ApJ, 738, 131 

Dutrey, A., Guilloteau, S., & Guelin, M. 1997, A&A, 317, L55 

Dutrey, A., Henning, T., Guilloteau, S., et al. 2007, A&A, 464, 615 

Dutrey, A., Wakelam, V., Boehler, Y., et al. 2011, A&A, 535, A104 

Duvert, G., Guilloteau, S., Menard, F., Simon, M., & Dutrey, A. 2000, A&A, 355, 
165 

Espaillat, G., Galvet, N., D’Alessio, P., et al. 2007, ApJ, 670, L135 
Fuente, A., Gernicharo, J., Agundez, M., et al. 2010, A&A, 524, A19 
Goldsmith, P. F., & hanger, W. D. 1999, ApJ, 517, 209 
Groesbeck, T. D., Phillips, T. G., & Blake, G. A. 1994, ApJS, 94, 147 



- 15 - 

Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385 

Henning, T., Semenov, D., Guilloteau, S., et al. 2010, ApJ, 714, 1511 

Isella, A., Perez, L. M., & Garpenter, J. M. 2012, ApJ, 747, 136 

Kastner, J. H., Hily-Blant, P., Rodriguez, D. R., Punzi, K., & Forveille, T. 2014, 
ApJ, 793, 55 

Kastner, J. H., Hily-Blant, P., Sacco, G. G., Forveille, T., & Zuckerman, B. 2010, 
ApJ, 723, L248 

Kastner, J. H., Qi, G., Gorti, U., et al. 2015, ApJ 

Kastner, J. H., Zuckerman, B., Hily-Blant, P., & Forveille, T. 2008, A&A, 492, 469 

Kastner, J. H., Zuckerman, B., Weintraub, D. A., & Forveille, T. 1997, Science, 277, 
67 

Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 
Kraus, A. L., & Ireland, M. J. 2012, ApJ, 745, 5 

Mangum, J. G., & Shirley, Y. L. 2015, ArXiv e-prints, arXiv:1501.01703 

Maret, S., Hily-Blant, P., Pety, J., Bardeau, S., & Reynier, E. 2011, A&A, 526, A47 

Milam, S. N., Savage, G., Brewster, M. A., Ziurys, L. M., & Wyckoff, S. 2005, ApJ, 
634, 1126 

Mulders, G. D., Dominik, G., & Min, M. 2010, A&A, 512, All 

Muller, H. S. P., Klaus, T., & Winnewisser, G. 2000, A&A, 357, L65 

Oberg, K. L, Qi, G., Fogel, J. K. J., et al. 2010, ApJ, 720, 480 

Pietu, V., Dutrey, A., & Guilloteau, S. 2007, A&A, 467, 163 

Pietu, V., Dutrey, A., Guilloteau, S., Ghapillon, E., & Pety, J. 2006, A&A, 460, L43 

Pulliam, R. L., Edwards, J. L., & Ziurys, L. M. 2011, ApJ, 743, 36 



Qi, C., Kessler, J. E., Koerner, D. W., Sargent, A. I., & Blake, G. A. 2003, ApJ, 
597, 986 

Salter, D. M., Hogerheijde, M. R., van der Burg, R. F. J., Kristensen, L. E., & 
Brinch, C. 2011, A&A, 536, A80 

Schoier, F. L., Jprgensen, J. K., van Dishoeck, E. F., & Blake, G. A. 2002, A&A, 
390, 1001 

Scott, P. G., Asplund, M., Grevesse, N., & Sauval, A. J. 2006, A&A, 456, 675 

Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034 

Skatrud, D. D., de Lucia, F. G., Blake, G. A., & Sastry, K. V. L. N. 1983, Journal 
of Molecular Spectroscopy, 99, 35 

Skinner, S. L., & Giidel, M. 2013, ApJ, 765, 3 

Thalmann, G., Mulders, G. D., Hodapp, K., et al. 2014, A&A, 566, A51 

Thi, W.-F., van Zadelhoff, G.-J., & van Dishoeck, E. F. 2004, A&A, 425, 955 

Thi, W. F., van Dishoeck, E. F., Blake, G. A., et al. 2001, ApJ, 561, 1074 

van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 330, 145 

van Zadelhoff, G.-J., van Dishoeck, E. F., Thi, W.-F., & Blake, G. A. 2001, A&A, 
377, 566 

Watanabe, Y., Sakai, N., Lindberg, J. E., et al. 2012, ApJ, 745, 126 

Williams, J. P., & Gieza, L. A. 2011, ARA&A, 49, 67 

Zhu, Z., Nelson, R. P., Hartmann, L., Espaillat, G., & Galvet, N. 2011, ApJ, 729, 47 

Ziurys, L. M., Saykally, R. J., Plambeck, R. L., & Erickson, N. R. 1982, ApJ, 254, 
94 


This preprint was prepared with the AAS lATRX macros v5.2. 



17 


Table 1. LkCa 15: Molecular Species Detected in the 
IRAM 30 M Line Survey 


Species 

Transition 

V 

(GHz) 

Tmf) 

(K) 

Vd “ 

(km s“^) 

I 

(K km s“^) ^ 


J = 2 ^ 1 

230.5380000 

0.542 

1.128 

1.441 

13 CO 

J = 2^ 1 

220.3986765 

0.201 

1.184 

0.550 

Ci ®0 

J = 2^ 1 

219.5603568 

0.038 

1.124 

0.092 

CN 

= 2 ^ 1 

226.6321901 

0.036 

1.413 

0.108 



226.6595584 

0.079 

1.084 

0.226 



226.6636928 

0.040 

1.228 

0.118 



226.6793114 

0.042 

1.251 

0.112 



226.8741908 

0.183 

1.431 

0.681'’ 



226.8747813 






226.8758960 






226.8874202 

0.020 

1.762 

0.081 



226.8921280 

0.015 

1.188 

0.062 

cs 

J = 5^4 

244.9356435 

0.064 

1.245 

0.182 

G2H 

N = 3^2 

262.0042266 

0.049 

1.128 

0.269^ 



262.0064034 






262.0648433 

0.033 

1.128 

0.152^ 



262.0673312 




HGN 

J = 3^2 

265.8861800 

0.238 

1.411 

0.732 

HCO+ 

J = 3 ^ 2 

267.5575260 

0.334 

1.298 

1.024 

DON 

J = 3^2 

217.2386307 

0.019 

1.128 

0.152® 

DGO+ 

J = 3 ^ 2 

216.1125766 

0.015 

1.128 

0.118® 

H2CO 

j = 3i3 — >■ I12 

211.2114680 

0.056 

1.300 

0.147 


j = 3o3 2o2 

218.2221920 

0.032 

1.134 

0.085 


J = 322 — >■ 221 

218.4756320 

<0.008 

1.128 

<0.060 


J = 3i2 — >■ 2ii 

225.6977750 

0.058 

1.093 

0.136 


a) As defined by Kastner et al. (2008). 

b) Sum of integrated intensities of hyperfine structure lines in range 226.87419-226.87590 GHz. 

c) Sum of integrated intensities of hyperfine structure lines in range 262.00426-262.00649 GHz. 

d) Sum of integrated intensities of hyperfine structure lines in range 262.06498-262.06747 GHz. 

e) Tentative detection. 

f) These line intensities can be converted to Jy km s“^ using 3.906^^ Jy/K for S/T^, the point 
source sensitivity, with Ff.ff= 0.92 and A^ff = 0.46 (aperture efficiency), where these represent 
mean values for the IRAM 30 m telescope in the 1.3 mm window. 
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Table 3. LkCa 15: Molecular Transitions Observed to Date 
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a) Beam-averaged column densities assuming optically thin emission and adopting source radii (for all molecules) of 4" See § 3.3 

b) column densities derived directly from the intensity assuming [^®0d®0]/[^^Cd^C] = 7 (Scott et al. 2006: Milam et al. 2005 . 
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Table 6. Column Density Comparisons 


Species 

This Study 
Nx 

(cm“^) 

T 

ex 

(K) 

Previous Studies 

Nx Tex 

(cm“^) (K) 

Reference “ 

13CO 

9.97 X 10^3 

30 

1.11 X IOI 6 

30 

2 


6.18 X 10^4 

20 

3.6 X 1014 

25 

3 


7.28 X 

10 

2.6 ±0.2 X 1013 

13.7 

4 

CiSQ 

9.97 X 10^3 

30 

3.31 X 1013 

30 

2 

cs 

1.20 X 10^2 

20 

5.1 X IO 12 

25 

3 


3.65 X 10^2 

10 

< 1.2 X 1013 

10 
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2.49 X 1013 
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8.7 ± 1.6 X IO 12 

8.7 
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HCN 

1.16 X IO 12 

20 

1.8 X IO 12 

25 

3 


2.18 X IO 12 

10 

4.3 X IO 12 

10 

6 


7.84 X IO 12 
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10.6 ± 1.5 X IO 12 

7.0 ±0.6 
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HCO+ 

9.12 X lOii 

30 

9.25 X IO 12 

30 
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9.47 X lOii 

20 

3.3 X lOii 

25 
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9.47 X lOii 

20 

8.0 ±0.5 X IO 12 
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1.80 X IO 12 
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DCO+ 

< 1.17 X lOii 
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< 2.9 X IQii 
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CN 

(0.52 - 0.55) X 1014 

3.9 ±0.2 
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25 
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(0.52 - 0.55) X 1014 

3.9 ±0.2 
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10 
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(0.52 - 0.55) X 1014 

3.9 ±0.2 

58±5 X IO 12 

8.8 ±0.3 
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C 2 H 

(1.29-2.01) X lOi® 

2.9 ±0.1 

2.9 ± 1.1 X 1013 

6.3 ± 1.4 

5 


a) References: 1 

. Aikawa et al. (20031; 2. Qi et al. 

(2003); 3. Thi et al. (2004); 4. Pietu et al. 

to 

0 

0 

); 5. 

Henning et al. 

(2010) 

6. 

Fuente et al. 

(2010); 7. Dutrey et al. 

(20111; 8. 

Chapillon et al. (2012a) 











































26 


O 

O 

CO 

§ 

s 


o' 

o 


o 

H 

w 

> 


w 

m 

W 

D 

p 

P 

m 

p 

< 

o 

P 

o 

< 

Pi 

[Pi 

p 

o 

z 

o 

S 

< 

p 

s 

o 

o 


1 ^ 


^ CD 
03 ^ 

5^ g 

o + 

■S o 

6 ^ 


CO 

a: 

CM 


O 

zn 

>- 


!X! 

Q 

i>^ 

a 

o 

S 

3h 

O 

o 

pH 


3 

CJ 

o 


I:-- 

co" 


zn ^ 

«< 

Pi 


m 

u 

Pi 


m 

M CO 

<; 


Pi 


< 

m CO 

<; 


Pi 


'oO 

m 

CD 

'Ct^ 

o 

'ct^ 

> 


cS 

M 


c3 

o 


CO 

kO 


kC 


zo 

zd 00 


1—I 1—I OP 

o o o 
o o o 


GO 00 CO 

1—I TJ* t—I 


'rt< Tt* CO 

( I I I I ( 

CO 00 ^ CO ‘O CN 

kO kO 00 OP OP CO 


CO 

'dy '~~L^ ■ '1' 

T—I CO T—i OP 
CO CO 


CO CO CO 

'dy '~~L^ ■ '-L' 

^ ^ CO ^ 
^ ^ 


CO 

I 

t-- 


co 

^ 00 
O ^ 00 

n: o 

V 


CO 

CO 

o 


sj^ 


CO i 

! 00 

o 

CM 

kO 



o 


G 



c? 



-S 

- S 


-V 


CO 

CO 

CO 

CO 

o 

1 

sj^ 


sj^ 


CO 

CM 

kO 

CO 

00 

kO 

OP 

p 

H 

T—1 

t-H 

CO 

CO 

o 

1—1 


CO 


+ O HH 

R ^ o M ^ 

p « o o o 


K M 





2; 

o 


^ CO 

o 


bO 


T3 

o 

o 

a 

s 

o 

CO 

CO 

o 

CO 

d 

X3 

-rs 

a 




CO 

OP 

(M 

CO 


zn 

< 


o ^ 

CO 

kC 

o ^ 

!> 

H 

.2 ^ 

bS) 


X3 

o 

a ^ 

I ^ 

c3 Ct3 
CO 

5 CO 

rO; ^ 

CM 

> 9 
































27 



Frequency (MHz) 


Fig. 1.— Receiver temperatures (top), system temperatures (middle), and signal- 
band opacities (bottom) as functions of frequency during the LkCa 15 spectral survey. 
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Fig. 2.— IRAM EMIR spectral coverage for the LkCa 15 spectral survey. 
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Fig. 3.— Spectra of all molecular transitions that have been studied in our IRAM 
30m line survey of LkCa 15. Ordinate is velocity with respect to the local standard 
of rest (in km s“^) and abscissa is line flux (in K); for clarity, spectral regions (other 
than ^^CO) have been shifted upwards, and lines have been rescaled. 
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Fig. 4.— IRAM 30m spectra of (J = 2 —)■ 1) (top), (J = 2 —)■ 1) (middle), 
and (J = 2 —)■ 1) (bottom) detected in LkCa 15 (blue histogram) overlaid with 
parametric Keplerian model line profiles (see § 3.2). Ordinate is velocity with respect 
to the local standard of rest (in km s“^, with frequencies in MHz at the top of the 
panel) and abscissa is line flux (in K). 


Frequency (MHz) 

230545 230540 230535 230530 


: 1 1 1 1 1 1 1 r 

CO (2-1) 

: J 


r : 

!1_1_1_1_1_1_1_1_1_L 



-5 0 5 10 15 

Velocity (km/s) 

Frequency (MHz) 


220405 220400 220395 220390 



Frequency (MHz) 

219565 219560 219555 






















































31 


Fig. 5.— IRAM 30m spectra of HCO’*' (J = 3 —)■ 2) (top), HCN (J = 3 —)■ 2) 
(middle), and CS (J = 5 —)■ 4) (bottom) detected in LkCa 15 (blue histogram) 
overlaid with parametric Keplerian model line prohles (see § 3.2). Ordinate is velocity 
with respect to the local standard of rest (in km s“^, with frequencies in MHz at the 
top of the panel) and abscissa is line flux (in K). 
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Fig. 6.— IRAM 30m spectra of H 2 CO (J = 3i3 —)■ 212 ) (top), H 2 CO (J = 3o3 —)■ 202 ) 
(middle), and H 2 CO (J = 3i2 —)■ 2ii) (bottom) detected in LkCa 15 (bine histogram) 
overlaid with parametric Keplerian model line prohles (see § 3.2). Ordinate is velocity 
with respect to the local standard of rest (in km s“^, with freqnencies in MHz at the 
top of the panel) and abscissa is line flnx (in K). 
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Fig. 7.— IRAM 30m smoothed spectra of DCN (J = 3 —)■ 2) (top) and DCO+ 
(J = 3 —)■ 2) (bottom) detected in LkCa 15 (bine histogram) (see § 3.2[ ). Ordinate 
is velocity with respect to the local standard of rest (in km s“^, with freqnencies in 
MHz at the top of the panel) and abscissa is line flnx (in K). 
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Fig. 8.— Full spectra of liyperfine components of CN observed toward LkCa 15. The 
red lines indicate the positions and relative intensities of the hyperhne components 
(See Table [^. The spectra are overlaid with the best-£t models obtained via the 




LSR velocity (km/s) 




















- 35 - 


Fig. 9.— Full spectrum of hyperfine components of C 2 H observed toward LkCa 
15. The red lines indicate the positions and relative intensities of the hyperhne 
components (See Table|^. The spectrum is overlaid with the best-£t models obtained 
via the method described in § |3.5[ The residuals from the fit are below the spectrum. 
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